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Conjugate  Silhouette  Nets 

W.  L.  F.  Degen 


Abstract.  Conjugate  nets,  Laplace  transformations  and  projective 
translation  surfaces  are  exploited  for  CAGD  purposes.  The  latter  are 
shown  to  be  equivalent  with  conjugate  nets  having  degenerated  Laplace 
transforms.  Relations  to  conjugate  nets  with  planar  silhouettes,  supercy- 
clides  and  Dupin  cyclides  are  given. 


§1.  Conjugate  Nets  and  their  Laplace  Transforms 

Conjugate  nets  play  an  important  role  in  classical  differential  geometry,  espe- 
cially because  of  their  projective  invariance  (see  [3]).  Representing  a surface 
in  d-space  by  homogeneous  coordinates 

<5  ...  nfl^R^flcR2,  (1) 

where  D is  an  open  connected  domain  of  the  “parameter  plane”  It2,  then  a 
conjugate  net  is  defined  by  the  validity  of  a Laplacian  equation 

xuv  + axu  + bxv  + cx  = 0 (2) 

with  certain  functions  a,b,c  € C°°\D\.  To  exclude  planar  surfaces 

dim sp&n(x,xu,xv,xuu,xvv)  ^3,  d ^ 3 (3] 

is  assumed  throughout.  The  geometric  meaning  of  (2)  for  d = 3 is  that 
the  characteristic  lines  of  the  tangent  planes  along  one  isoparametric  line  are 
tangent  to  the  other.  A second  characterization  is  that  the  two  tangents  to 
the  isoparametric  lines  have  a harmonic  cross  ratio  with  the  two  asymptotic 
tangents. 

The  most  important  property  of  conjugate  nets  is  that  they  have  in  both 
isoparametric  directions  a Laplace  transformed  net  ( Laplace  transform)  which 
is  defined  by  the  striction  points  of  the  two  developables  generated  by  the 
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tangent  planes  along  the  net  curves  (isoparametric  lines).  Analytically,  they 
are  given  by 

£„(*)  •••  y:=xu  + bx,  (4) 

Cv(x)  ...  z:=xv  + ax.  (5) 

Indeed,  by  (2),  we  get 

yv  = hx-  ay, 
zu  = kx  — bz 

with 

h = bv  4-  ab  — c,  k = au  + ab  — c (8) 

being  the  so-called  Darboux  invariants.  Eqn.  (6)  shows  that  the  line  x A y is 
tangent  to  the  u-curve  (isoparametric  line  with  v = const,  u varying)  on  the 
surface  (x) D as  well  as  tangent  to  the  u-curve  on  (y) p ; analogously  x A z is 
tangent  to  the  u-curve  on  (x)p  and  the  u-curve  on  (z)D. 

In  euclidean  differential  geometry,  conjugate  nets  owe  their  importance 
to  the  fact  that  the  curvature  lines  on  every  surface  are  conjugate  (except 
at  umbilic  points).  Thus,  for  CAGD  applications,  one  can  profit  from  the 
theory  of  conjugate  nets  since,  for  instance,  the  two  families  of  circles  on  a 
Dupin  cyclide  or  the  net  of  profile  curves  and  meridian  circles  on  a surface  of 
revolution  is  a conjugate  net. 

In  this  paper  we  deal  with  the  special  case  that  the  two  Laplace  transforms 
degenerate  into  curves.  Assuming  that  the  net  curves  have  nowhere  asymptotic 
directions,  i.e. 


(6) 

(7) 


x A xu  A xuu  ^ 0 , x A xv  A xvv  ^ 0,  (9) 

then  the  degeneration  conditions  y A yu  A yv  = 0 and  z A zu  A z„  =0  imply 

h = 0,  k = 0 for  all  (u,v)  £ D.  (10) 

Definition  1.  A conjugate  net  in  3-space  satisfying  (10)  with  regularity  con- 
ditions (3),  (9)  is  called  a conjugate  silhouette  net. 

This  notation  is  justified  since  by  (4),  (6),  (10)  all  w-tangents  along  a 
u-curve  meet  at  the  fixed  point  y,  thus  building  up  a general  cone  with  apex 
y.  Thus  the  u-curve  x(u0)  v)  can  be  considered  as  a silhouette  on  the  surface 
S by  central  illumination  from  y.  Similarily  the  u-curves  are  silhouette  lines 
by  central  illumination  from  z. 

On  the  other  hand,  it  is  easy  to  see  that  a net  of  silhouette  lines  on  a 
surface  whereby  the  centers  of  illumination  vary  on  two  curves  is  automatically 
conjugate,  provided  that  the  generators  of  the  enveloping  cones  are  tangent  to 


Conjugate  Silhouette  Nets 


39 


the  net  curves.  Thus,  the  equations  (2),  (10)  characterize  conjugate  silhouette 
nets  (up  to  degenerated  cases). 

At  first  glance,  this  class  of  surfaces  seems  to  be  very  restricted.  But 
this  is  not  true:  it  comprises  many  subclasses  of  surfaces  considered  in  CAGD 
literature  such  as  Dupin  cyclides  ([2,4,7]),  supercyclides  ([1,5, 8, 9]),  nets  with 
planar  silhouettes  ([6])  etc..  As  for  all  of  these  examples,  the  whole  class  of 
conjugate  silhouette  nets  is  well-suited  for  CAGD  applications,  in  particular 
for  geometric  modelling  purposes  because  of  their  simple  blending  properties: 
Putting  two  of  them  together  along  a common  net  curve  immediately  yields 
a (^-continuity,  once  the  corresponding  centers  of  illumination  coincide. 

But  there  is  still  another  reason  making  these  surfaces  worth  consider- 
ing in  CAGD:  They  admit  a very  simple  generation  as  so-called  “projective 
translation  surfaces”,  as  will  be  derived  in  the  next  section. 

§2.  Projective  Translation  Surfaces 
Let 


A ...  p : Ii  —y  lRd+1,  £2  ...  q:I2-+Rd+1  (11) 

be  two  C^-curves  in  d-space  represented  also  in  homogeneous  coordinates 
(Ii,  I2  being  two  open  nonvoid  intervals  of  1R).  Then  one  gets  a surface  S (1) 
simply  by  setting 

S . . . x(u,v)  :=  p(u)  + q(v),  (u,v)  E Ix  X I2  D (12) 

Definition  2.  Surfaces  defined  by  (12)  via  two  curves  (11)  are  called  projective 
translation  surfaces. 

This  definition  generalizes  the  usual  euclidean  (or  affine)  definition  of 
translation  surfaces,  where  the  same  formula  (12)  is  used  but  interpreted  in 
affine  (non-homogeneous)  coordinates.  So  one  curve  can  be  considered  to 
move  along  the  other  thus  sweeping  out  the  surface.  In  the  projective  case, 
the  generating  point  x(u,  v ) always  lies  on  the  line  p(u)/\q(v)  joining  these  two 
points  of  Ci  and  C2  independently.  It  must  be  noticed  that  the  normalizations 
are  essential  (not  arbitrarily  to  be  choosen  like  usually  when  dealing  with 
curves):  they  determine  the  position  of  that  point  x(u,v)  on  the  line  p(u)  A 
q(v). 

Now  we  can  establish  one  of  our  main  results: 

Theorem  1.  Every  conjugate  silhouette  net  is  a projective  translation  sur- 
face, and  the  net  curves  correspond  to  the  isoparameter  lines  in  the  represen- 
tation (12). 

Proof:  We  have,  by  definition,  h = 0,  k = 0,  and  hence  in  particular 


a y — by  . 


(13) 
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Assuming  D to  be  simply  connected,  we  conclude  that  there  exists  a C°°- 
function  /:£>—>  1R  with 

fu  =b , /„  = a.  (14) 

Taking  p :=  e? , we  calculate  pu  = pb,  pv  = pa  and  puv  = ( pb)v  = p(ab  + bv)  = 
pc,  the  latter  observing  (8),  (10).  Since  p ^ 0 in  D,  we  can  renormalize 
x :=  px,  obtaining 

Xnv  = (Puv  - pc)x  + ( pu  - pb)xv  + (pv  - pa)xu , (15) 

and  thus 

x„„  = 0.  (16) 

This  equation  immediately  yields  a representation  (12)  by  integration  (possi- 
bly restricted  to  a rectangle  I\  x I2  within  D).  □ 

Up  to  now  we  think  of  that  renormalization  as  always  having  been  done, 
so  the  Laplace  equation  (2)  has  the  coefficients 

o = 0,  5 = 0,  c = 0.  (17) 

Therefore,  the  Laplace  transforms  (4)  and  (5)  are  now  given  by 

Cu(x)  ...  *(u)  = ^*r.  »€/!,  (18) 

£»(*)  •••  z(v)  = — v € h-  (19) 

Calling  these  curves  the  projective  hodographs  of  p and  q respectively,  we 
can  state 

Corollary.  The  Laplace  transforms  of  a conjugate  silhouette  net  C are  the 
projective  hodographs  of  the  generating  curves  C\ , C2  of  C (considered  as  a 
projective  translation  surface). 

§3.  Axial  Silhouette  Nets 

Definition  3.  A conjugate  silhouette  net  is  called  axial  if  the  generating 
curves  C 1,  C2  in  its  representation  (12)  as  a projective  translation  surface  are 
(parts  of)  straight  lines.  These  lines  are  called  the  first  and  the  second  axis  of 
the  net. 

The  conditions  for  axial  conjugate  silhouette  nets  are  that  p,p',p'!  and 
likewise  q,  q' , q " must  be  linearly  dependent  (a  prime  at  p indicating  deriva- 
tion with  respect  to  u and  at  q with  respect  to  v).  Assuming  p,  p'  and  likewise 
q,  q'  to  be  linearly  independent  (otherwise  the  point  would  be  stationary)  we 
have 

p"  = ap  + Pp',  q"  = 7 q + Sq'  (20) 

with  some  C'°°-functions  a,/?  of  u and  7, 6 of  v characterizing  axial  nets. 

This  has  many  consequences;  most  of  them  we  proved  earlier  for  super- 
cyclides  and  for  nets  with  planar  silhouettes  [5,6].  Now  we  give  the  result  for 
the  general  case  of  axial  silhouette  nets: 
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Theorem  2. 

a)  All  the  net  curves  (of  both  families)  are  planar  curves, 

b)  The  planes  of  the  net  curves  of  each  family  belong  to  a pencil, 

c)  The  axes  of  these  two  pencils  coincide  with  the  second  and  the  Grst  axis 
of  the  net  (i.  e.  the  plane  of  a curve  of  the  Grst  [seconcf]  family  passes 
through  the  second  [first]  axis, 

d)  The  apexes  of  the  envelopping  cones  ( the  “light  centers”)  along  a u-curve 
[v-curve]  lie  on  the  second  [first]  axis, 

e)  Any  two  u-curves  [v-curves]  are  projectively  equivalent  to  each  other. 

Proof:  We  perform  the  proofs  only  for  the  u-curves;  the  assertions  with 
respect  to  the  e-curves  follow  analogously. 

a) :  From  (12),  (20)  we  derive 

xu  = P',  xuu  =ap  + (Ip1,  (21) 

xUuu  = (<*'  + a0)p  + (a  + /3'  + 02)p' . (22) 

Hence  xu  A xuu  A xuuu  = 0,  meaning  that  the  u-curves  are  planar. 

b) ,  c):  Eqns.  (21)  show  that  p,  p'  are  contained  in  the  plane  x AxuA  xuu 
of  this  u-curve;  but  p,  p'  span  the  first  axis. 

d) :  The  apex  of  the  envelopping  cone  is  given  by  (21)  as  p'.  Hence  it  is 
lying  on  the  first  axis. 

e) :  Assumption  (9)  implies  a ^ 0 for  all  u € h-  Thus,  xu  and  xvu  can 
be  eliminated  from  (21)  and  with  this  Eqn.  (22)  yields 

Xuu  d" 

Thus  the  coefficient  of  the  fundamental  equation  (see  [3])  do  not  depend  on  the 
second  parameter  v;  this  means  geometrically  that  all  u-curves  are  projectively 
equivalent.  □ 

For  CAGD  purposes  rational  (and  polynomial)  surfaces  are  of  particular 
interest.  The  explicit  representation  (12)  makes  it  very  easy  to  pick  out  ratio- 
nal surfaces  from  that  general  class:  The  only  thing  one  has  to  do  is  to  insert 
rational  representations  for  p(u)  and  q{v).  We  restrict  this  procedure  to  axial 
nets,  and  derive  from  it  the  (rational)  Bezier  representation. 

Theorem  3.  For  any  pair  of  planar  rational  curves  C\  and  C2,  there  exists 
an  axial  conjugate  silhouette  net  having  its  u-curves  projectively  equivalent  to 
C 1 and  its  v-curves  projectively  equivalent  to  C 2.  The  axes  can  be  arbitrarily 
prescribed  as  two  skew  straight  lines. 

Proof:  Let  Ci,  C 2 be  represented  in  planar  homogeneous  coordinates  by 
triples  of  linearly  independent  polynomials 


a 
, 0 


0 + a + 0' ) xu. 


(23) 


IUU  — [ P 

a 


C. 1 ...  Xi  = fi(u)  ( i = 0, 1, 2),  C2  . . . Xi  = gi(v)  (i  = 0, 1, 2).  (24) 
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Furthermore,  let  the  axes  be  spanned  by  vectors  a,  6 € ft4  and  c,  d € R4 
resp.  Then  det(a,  b,  c,d)  ^ 0 since  the  axes  are  assumed  to  be  skew.  With 
this  we  can  set 

P{u)  ~ j^{h{u)a  + f2{u)b),  q(v)  :=  -~^(gi(v)c  + g2{v)d)  (25) 

getting  the  desired  axial  net  S by  (12)  (restricted  to  intervals  /i,/2  C R 
where  fo(u)  resp.  go(v)  have  no  zeros).  Now,  indeed,  the  w-curves  are  planar 
and  projectively  equivalent  to  C\  since,  for  fixed  v = vo,  we  have  p(u)  = 
fo(u)p(u)  = fo(u)q(vo)  + fi(u)a  + f2(u)b  so  that  fi(u)  are  the  coordinates 
with  respect  to  basis  q(v o),  a,  b.  The  proof  for  the  u-curves  follows  the  same 
line,  mutatis  mutandis.  □ 

Obviously,  the  representation  (25)  is  not  unique.  But  we  can  immediately 
derive  from  (12)  and  (25)  the  Bezier  representation  of  S:  First  renormalize 
(12)  with  the  factor  fo{u)go(v)  getting 

x(u,v)  = g0(v)(fi(u)a  + f2{u)b)  + f0(u)(gi(v)c  + g2{v)d) , (26) 

and  then  we  expand  the  polynomials  fi(u)  and  Pi(v)  with  repsect  to  the  Bern- 
stein basis 


/<(«)  = 9k(v)  = '£^iBr(v),  (27) 

j=0  1=0 

getting  the  usual  homogeneous  Bezier  representation 

n m 

*(«.«)  = B?(v)  (28) 

j= 0 !=0 

with  the  control  points 

bj,i  = /30,i(«i,jO  + «2  jb)  + aojiPijc  + 02,id).  (29) 

Since  these  calculations  can  be  done  also  backwards,  we  obtain 

Corollary.  The  conditions  (29)  for  the  control  points  characterize  a (n,  m )- 
rational  Bezier  surface  (28)  to  be  an  axial  conjugate  silhouette  net. 


§4.  Applications  to  Dupin  Cyclides 

Dupin  cyclides  are  special  kinds  of  supercyclides  [5],  and  therefore  they  should 
have  representations  as  axial  conjugate  silhouette  nets  so  far  they  have  fourth 
order  and  skew  axes.  However,  their  usual  representation  starts  from  the 
midpoint  curves  Y(u)  and  Z(v)  of  the  two  families  of  envelopping  spheres 
and  results  in 


V ...  X(u,v)  = 


r2(v)Y(u)  + ri(u)Z(v) 
r2(v)  + rj(u) 


(30) 
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where  ri,r2  denote  the  corresponding  radius  functions.  Though  this  is  also  a 
sum  of  two  vectors  it  has  neither  the  form  (12)  nor  are  the  components  parts 
of  straight  lines. 

Thus,  the  question  arises  of  how  to  transform  (30)  into  (12).  The  first 
step  to  solve  this  problem  consists  in  passing  to  homogeneous  coordinates 
x = (x0,xi,x2,x3)t 

* = pQ-),  peiR\{o}  (31) 


(i.e.  xq  = p • 1,  X{  = p ■ Xi  (i  = 1,2,3))  and  to  take  in  the  present  case 
P = — + ~ yielding  indeed  (12) 


/ ) ' \ 

( 1 \ 

n(u) 

, q(v)  = 

r2(v) 

l 

1 

\r2(u)Z{u)J 

(32) 


However,  the  curves  p and  q describe  an  ellipse  and  a hyperbola  as  before. 
Obviously,  another  representation  of  the  same  kind  (12)  must  satisfy 


p{u)  = p(u)  + c,  q(v)  = q{v)  - c (33) 

with  a constant  vector  c. 

Starting  with  the  explicit  representations 


/ 1 \ 

( 1 ^ 

1-u2 

1 + V2 

Y(u)  = 

l + «2p 

, Z(v)  = 

i 2 P° 

1 — V1 

2u 

0 

1 + u*PUJ 

2v 

V o y 

(34) 


(with  some  shape  parameters  d,p,o,  [cr|  < 1,  w = \/l  — c2)  and  observing  the 
corresponding  radius  functions 


n(u) 


fo{u)  gojv) 

1 + u2  ’ 2 1 — V2  ’ 


(35) 


whereby 


fo(u)  - (1  + u2)dp  - (1  - u2)pcr,  g0(v)  = (1  + v2)p  - (1  - v2)do, 
we  finally  obtain 

/ /1\  ( 0 \\ 


p(u) 


fo(u) 


((!  + u2)  - — f0(u)) 


\ 


d 

0 

voy 


+ 


0 

2 puu 

\ o y 


(36) 
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and 


q{u) 


9o(u) 


((i  - v2)  + — g0(v)) 


V 


/ 1\ 
da2 
0 

V 0/ 


+ 


/ 0 \ 
0 
0 

\ 2puiv  J 


(37) 


with 


Thus  we  proved 


e=£( i»  o.  °>  °)T- 


Theorem  4.  The  formulas  (36),  (37)  (inserted  into  (12))  yield  an  explicit 
representation  of  nonparabolic  Dupin  cyclides  with  skew  axes  as  axial  con- 
jugate silhouette  nets  with  respect  to  a suitable  (homogeneous)  coordinate 
system. 
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